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We consider a plasma described by means of a two-dimensional fluid model across a constant but 
non-uniform magnetic field B = B(x,y)z. The dynamical evolution of the density and the vorticity 
takes into account the interchange instability and magnetic field inhomogeneities. First, in order 
to described the Finite Larmor Radius effects we apply the gyromap to build a Hamiltonian model 
with ion temperature from a cold-ion model. Second, we show that the gyromap is justified using 
Braginskii's closure for the stress tensor as well as an apt ordering on the fluctuating quantities. 

I. INTRODUCTION 

In laboratory plasma devices like tokamaks or in astrophysical plasmas, two-dimensional reduced fluid models have 
been used to investigate various instabilities and their impact on turbulent transport. These reduced models are very 
useful in practice since their theoretical and numerical analyses are simpler than for complete models. At the same 
time they display important features of the underlying turbulent phenomena with more clarity than the complete 
models. The standard procedure used for the derivation of these reduced models is to start from the equations of 
motion of a parent model (e.g., that given by the continuity and the momentum equations) and to perform an ordering 
with respect to one or several small parameters on the variables of the system (e.g., density and potential fluctuations) 
and on the electromagnetic field. This ordering is suggested by the physics and the geometry of the phenomenon under 
consideration (e.g., one might want to filter out irrelevant time or length scales). It is desirable that the reduced model 
preserves the main ingredients of the parent model under this reduction procedure, as for instance the Hamiltonian 
symmetry. Typical parent models have indeed this Hamiltonian symmetry, e.g., the collisionless Maxwell- Vlasov 
equations [3, 0] or ideal magnetohydrodynamics However, in such a derivationprocedure obtained by working 
at the level of the equations of motion, neglecting terms can produce "mutilations" [4j of the original fluid equations 
by breaking this symmetry. These mutilations, which introduce incorrect dissipative terms, have a drastic impact 
on the properties of the physical system. The loss of the Hamiltonian structure indeed generates a more intricate 
interpretation of the numerical results, and in general clouds the relation with the parent model and its properties. 
In fact, as shown for example in models describing magnetic reconnection 0, Q, the Hamiltonian structure helps the 
interpretation and explanation of complex physical processes observed in numerical simulations. Among the other 
advantages of the Hamiltonian structure, there is the possibility of identifying systematically invariants of motion 
(e.g., Casimir invariants or conserved quantities linked to continuous symmetries), and applying tools of Hamiltonian 
perturbation theory and methods to investigate the stability of equilibria. 

However, for some reduced models, the important property of possessing a Hamiltonian structure has not been shown. 
For instance, this is the case of the two-dimensional fluid model taken from Ref. Q ■ This four-field model includes 
Finite Larmor Radius (FLR) effects as in Refs. 0,0], the drift velocity ordering as in Ref. [l(| and the gyroviscous 
terms as in Refs. [HI, Gil- Moreover, magnetic inhomogeneities are kept in the continuity equation as in Refs. fj\. [l3l— 
Il5l | . For this particular model, the Hamiltonian structure has not been found even though the model has been shown 
to be energy conserving. In this paper, we consider the main obstacle associated with the search for the Hamiltonian 
structure of this four-field model which is already encapsulated in a two-field model for the density and the vorticity 
fields. This two-field model is obtained from the four-field model by suppressing the parallel dynamics (along the 
magnetic field), and the poloidal magnetic and parallel flow fluctuations, and reads 



On 
dt 



1 

B 



b,A<f,] + (l + T) 



1 

B 



n ' ~B 



TV • [n, V<j 



(1) 



(2) 



where n = n(x, y, t) is the logarithm of the normalized density fluctuations, <f> = <t>(x, y, t) is the normalized electrostatic 
potential, B — B(x, y) is the normalized magnetic field, T is the constant ion temperature normalized by the electron 
temperature and [/, g] = z ■ V/ x Vg is the canonical bracket in the plane across the magnetic field B = Bz. All the 
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derivatives arc defined on the perpendicular plane. The bracket [<f>, n] corresponds to the dynamics due to the E x B 
drift, and the bracket [n, 1/B] corresponds to the term driving the interchange instability. 

When addressing the question of the Hamiltonian structure of the model ([I])-©, it is easy to show that, provided 
that boundary terms vanish, a conserved quantity for the model is 
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where the first part corresponds to the potential energy for ions and electrons and the second part is the kinetic 
energy. This conserved quantity / is interpreted as the total energy of the system and is a natural candidate for the 
Hamiltonian. Nevertheless, in spite of having a conserved quantity, the model (HJ-® docs not seem to possess a 
Hamiltonian formulation, at least of the Lie-Poisson type, which is the most common form for reduced fluid models. 
More precisely, Poisson brackets between functionals F and G of the form 
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where the vector bracket is defined by [A; B] = ^ [Ai, £?,], the terms a*, and are constant, £i = n , £2 = 
and SF/St; is the functional derivative of the functional F with respect to £, fail to give the desired equations of 
motion (d])-© when applied to /. Even if a bracket giving Eqs. ([I])-© with the invariant I defined by Eq. (J3|) has 
been found, it does not satisfy Jacobi identity and thus is not a Poisson bracket. 

The reason for the difficulty to find a Hamiltonian structure resides in the co-existence of compressibility terms (those 
containing 1/B in the continuity equation of the electrons) with the ion-gyroviscous term (in the vorticity equation). 
The compressibility terms are often retained in the continuity equation (see, e.g., Refs. 0,0413) i n order to account 
for the fact that the velocity field is not incompressible in the presence of an inhomogeneous magnetic field. If these 
compressibility terms are neglected in the two-field model, the resulting system is Hamiltonian (la , fToj ] . On the other 
hand, it has been shown that the system obtained by keeping these terms and eliminating the ion-gyroviscous terms 
is also Hamiltonian [l3| . It is thus the simultaneous presence of the two contributions which seems to complicate the 
search for a Hamiltonian structure. This leads to the problem of finding a Hamiltonian model accounting for both 
compressibility terms and FLR corrections. 

An elegant and effective way to introduce FLR corrections to a cold-ion model, while preserving the Hamiltonian 
structure, was introduced in Refs. [H|,[l7j and was referred to as the gyromap. The gyromap procedure was rigorously 
found from the Braginskii's closure for the stress tensor [l8| only for the three-field model described in Refs. [H, [l3[ 
which does not contain the compressibility terms in the continuity equation. In this paper, we apply the gyromap 
procedure to the cold-ion limit of the two-field model with the compressibility terms. We show that all the FLR 
correction terms produced by this method are obtained from the Braginskii's closure for the stress tensor, by means 
of an expansion based on a physically sound ordering. 

The paper is organized as follows: In Sec. 2, we apply the gyromap procedure to the cold-ion version of the two-field 
model ((I])-©. In Sec. 3, we show that the resulting model is directly obtained from a systematic expansion. 



II. THE GYROMAP 

The gyromap algorithm was introduced in Ref. (l3l | in order to introduce FLR corrections to a cold-ion model, while 
preserving the Hamiltonian structure. The model with FLR corrections obtained through this procedure possesses 
the same bracket of the original cold-ion model, but different dynamical variables and Hamiltonian. The fact that the 
Poisson bracket is inherited from the cold-ion Hamiltonian model ensures that the resulting model is Hamiltonian. 
As a by-product, the Casimir invariants of the resulting model with FLR corrections become easily available. We also 
recall that, as mentioned in Ref. [TH, the gyromap procedure possesses the desirable features that the cold-ion limit 
(T = 0) of the post-gyromap model gives the initial cold-ion model, and that the diamagnetic effects predicted by the 
kinetic theory are conserved at the first order. In what follows we use dimensionlcss quantities, i.e. we rescale space 
variables by the sonic Larmor radius p s , the density by the equilibrium density no, the charge by the electric charge 
e, the electron and ion temperatures by the electron temperature T e , the magnetic field B by the spatial mean value 
Bq and time by the inverse of the ion cyclotron frequency u) Cj i = e-Bo/m; where m\ is the ion mass. 
We start with the cold-ion version of Eqs. ((I])-© which describes the dynamical evolution of the plasma density n 
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and the vorticity A</> 
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We consider the system whose field variables are n(x,y) and 4>(x,y). Here we do not specify the time dependence 
of the variables which is implicitly assumed. In the algebra of observables which are functionals of n and </>, the 
Hamiltonian structure is defined by the Hamiltonian 
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and by the Poisson bracket 
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This bracket verifies all the required properties for a Poisson bracket like bilincarity, antisymmetry, Jacobi identity 
and Leibniz identity. 

The gyromap procedure is initiated from the Hamiltonian defined by 



H(n, $) = i J d 2 x ^(1 + T)n 2 + | V$| 2 ^j , 
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where <£> is the stream function for the FLR model. The choice of Eq. (jl]) as Hamiltonian is motivated by the 
requirement of adding the ion internal energy to the cold-ion Hamiltonian and having a kinetic energy term whose 
relation with the cold-ion kinetic energy is retrieved a posteriori. We introduce an auxiliary variable £ defined by 
A$ = £ + TAn/2. The shift of TAn/2 corresponds to a shift of half the magnetization velocity. This is the required 
transformation in order to yield the proper FLR corrections [l3|. The next step is to use the same Poisson bracket 
written in the new variables (iV, £), 
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The change of variables, N = n and £ = A$ — TAn/2, includes a change of functional derivatives, 
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Hence the bracket expressed in the variables (n, $), which is still a Poisson bracket, is the following, 
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With the Hamiltonian defined by Eq. Q, the equation of motion for the density is 

1 
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We assume that the continuity equation for electrons is unaffected by FLR corrections. Therefore comparing Eq. JT]) 
with Eq. © gives the following relation between $ and <f>: 
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On the other hand, the equation for the generalized vorticity A$ is given by 
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Equations ©-([9]) are the final result of the gyromap procedure. They include FLR corrections and the resulting 
system has a Hamiltonian structure since the bracket © is obtained from the Poisson bracket (J5J) by a change of 
variables. 

By expanding the operator (1 + TA/2) -1 for small T and making use of Eq. ([7]), it is possible to obtain an expression 
for Eq. ([9]) valid for small ion temperature T (in terms of n and </>), 
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This expression is useful, because it allows one to directly compare the leading order terms of the vorticity equa- 
tion ([TU| . with the vorticity equation as given in Ref. [l6| [see also Eq. ©] which refers to a model where the only 
FLR term is — TV - [n, V0]. We notice that already at the order O(T), the two equations differ from an FLR correction 
term given by — T [V</>; VA</>]. From Eq. (JTUJ), it is also straightforward to see that, in the limit T — > 0, the FLR 
model is correctly reduces to the original cold-ion model, as expected. 

In summary, in order to obtain a Hamiltonian structure with the ion temperature, the gyromap procedure starts from 
the Hamiltonian structure without ion temperature and then considers FLR effects using a change of variables. 
In what follows, we consider Braginskii's closure for the stress tensor in order to derive a model with FLR terms and 
compare it with the above model given by Eqs. (JTJ-©. 



III. DERIVATION FROM BRAGINSKII'S CLOSURE 



We first consider the dynamics of a plasma composed of two species, ions and electrons. In a simplified magnetic 
geometry where the magnetic field is constant (but non-uniform) and its direction is fixed, we restrict the plasma 
dynamics to the plane transverse to the magnetic field lines, i.e. all dynamical variables only depend on the two 
coordinates x and y. The dynamical variables are the ion and electron densities, n;(x, y) and n e (x, y), and the ion and 
electron velocity fields ~V{(x,y) and v e (x, y). We assume a quasi-neutrality hypothesis which allows us to reduce the 
number of variables, i.e. we consider that i%i(x, y) = n c (x, y) which is denoted n(x, y) in what follows. The dynamical 
equation for the density is given by the continuity equation for the electrons, 

| + V.(nv.) = 0. 

Due to the quasi-neutrality assumption, the same equation holds for ions, which translates into 

V ■ (n(v e - Vi)) = 0. (11) 
The dynamical equation for the velocity field v s (x,y) of the species s, where s refers to ions or electrons, is given by 

Hen + v s ■ V J v s = n (E + v s x B) - e s T s Vn - e s V • fj, (12) 



where E is the electric field, T s is the dimensionless temperature of the species s (i.e. T = T and T c = 1), /i s its 
dimensionless mass (i.e. /i; = 1 and u c = m c jm\ <g; /Zj because the mass of the electron m c is negligible compared to 
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the one of the ions mi) and e s its dimensionless charge (i.e. e\ — \ and e c = —1). The stress tensor associated with the 
species s, denoted 7fJ, is taken from Ref. [l8|. For each species, it is composed of viscosity terms identified by viscosity 
coefficients of various orders. Within a strong magnetic field approximation, we only consider the two higher order 
viscosity terms labeled by tjq and 773 with 770 3> 773. The viscosity coefficients for ions are t/q = 0.96nTri 3> 7? 3 = nT/2 
where t s is the normalized collisional time for the species s. 

In what follows, we only use the part of the stress tensor perpendicular to the magnetic field lines. Given these 
approximations, the divergence of the stress tensor for ions becomes 

_ T 
V • 7fT = -oiiTiTV (nV • vi) + — [nz x Av ; - (z X Vn • V)v ; + (Vn • V)(z x Vj)] , 

where ai = 0.32. The viscosity coefficients for electrons are t]q = 0.73nr c 3> ry| = — n/j, e /2 where r e is the collisional 
time for the electrons. We consider the following approximation for the divergence of the stress tensor for the electrons 

V • IzZ = -a c r c V (nV • v e ) , 
where a Q = 0.24. The dimensionless equation for the velocity field of each species becomes 

fi s n ( + v s • V J v., = -nV0 + nBv s X z - e s T s Vn + a s T s e s T s V (nV • v s ) 
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We apply the cross product with z/nB to Eq. (Ti~3"|) in order to obtain the velocity as the sum of the E x B, the 
diamagnetic, the polarization drifts and some gyroviscous contributions from the stress tensor, 
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We expand all dimensionless quantities in the following way 

l/B = l+e/B, 

(f) = £(j), 

11 = 1 + en, 

The small parameter e 1 denotes the amplitude of the fluctuations of the electromagnetic field as well as that of the 
dynamical variables n and v s . Moreover, the temporal variations of these quantities are small so that the temporal 
variation is denoted ed/dt. In what follows, we omit the tildas over the fluctuating quantities for simplicity. Here we 
consider the model for weak collisionality (i.e. /i s /r s <C 1) such that t; is of order 1/e and r c -C 1/e. In addition T 
is of order one, and the parameters a s are of order e (i.e. ol\t{F is of order 1 and a r c <c 1). The leading order of 
Eq. (TT4"]) (terms proportional to e) gives 

1 + Vs^^J vi 1 ' = z x V (0 + e s T s n) - a s T s e s T s z x V(V • v^ 1 '). 

From the previous equation, we directly obtain that V • v^ 1 ' = 0. From the expansion of Eq. (fTTj) . we also conclude 
that V • Vj = 0. Therefore the equation for Vg becomes 

l + Vs^Aj v« =zx V(0 + e s T s n), (15) 

so that the first order velocity is the sum of the E x B velocity and the diamagnetic velocity with FLR corrections 
on the left hand side. The second order of the velocity v s is given by 

1 + Ms^a) vi 2 ' = iv« + fj, s (JL + v« • v) z x vi 1 ) - eT s nz x Vn 

~a s T s e s T s z x V (V • v( 2 )) - [( z x Vn ■ V)(z x v^ 1 ') + (Vn ■ V)v( 1 )] , (16) 
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where we have used Eq. (TT5]) for Vg 1 '. We notice that we have also used the expansion V log( 1 + en) — e Vn — e 2 nVti + 
<3(e 3 ). 

As a first step, for the electrons, we see that from Eq. (|T5j) . the first order of the electron velocity is the sum of the 
E x B velocity and the electron diamagnetic velocity 
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and the second order of the electron velocity is 
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From this equation, we get that the divergence of is equal to the divergence of v^/B, and using Eq. (|TT)) . it 
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The continuity equation of the electron density at the second order in e becomes 
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or equivalently, using Eqs. (|T71) - (fT51) . 
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In fact, the first order of the ion velocity given by Eq. (|15|) suggests a change of variables (1 + -jA)$ 
which we can write the first order of the ion velocity into 



Tn with 



v> ' = z x V$. 

We notice that this change of variables is also the one suggested by the gyromap [see Eq. ©]. Using this variable $, 
the continuity equation becomes 
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As a second step, we work with the continuity equation for ions which is 

?P + V.v i (2) +v i (1) .V77 = 0. 
at 

First we obtain the following formula for the divergence of Vj from Eq. (1161) . 
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Next we multiply Eq. by (1 + TA/2) in order to use the Eq. (f!?Tj) and we insert Eq. P^|) so as to obtain 
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We notice that Eqs. (pT9 ]) -(|22 ]l coincide with Eqs. which were obtained by applying the gyromap. Therefore, 

we have shown that the terms generated from the gyromap are obtained consistently from Braginskii's closure for 
the stress tensor by making use of an appropriate ordering. 



In summary, it is possible to construct a model with FLR corrections and its Hamiltonian structure from a 
cold-ion model which possesses a Hamiltonian structure by applying a gyromap procedure which generates all the 
relevant FLR terms at the leading order. The change of variables introduced by the gyromap is directly given by 
the definition of the stream function at the first order of the ion velocity. We have shown that the two-field reduced 
model Q-© is obtained using the Braginskii's closure for the stress tensor by considering an apt ordering on the 
dynamical variables. 



7 



Acknowledgements 

We acknowledge financial support from the Agence Nationale de la Recherche. This work was supported by 
the European Community under the contract of Association between EURATOM, CEA and the French Research 
Federation for fusion studies. The views and opinions expressed herein do not necessarily reflect those of the European 
Commission. 



[1] P. J. Morrison, Phys. Lett. A 80, 383 (1980). 

[2] J. E. Marsden and A. Weinstein, Physica D 4, 058102 (1982). 

[3] P. J. Morrison and J. M. Greene, Phys. Rev. Lett. 45, 790 (1980). 

[4] P. J. Morrison, Phys. Plasmas 12, 058102 (2005). 

[5] D. Grasso, F. Califano, F. Pegoraro and F. Porcelli, Phys. Rev. Lett. 86, 5051 (2001). 
[6] E. Tassi, P. J. Morrison, D. Grasso and F. Pegoraro, Nucl. Fusion 50, 034007 (2010). 

[7] R. D. Hazeltine, M. Kotschenreuther and P. J. Morrison, Phys. Fluids 28, 2466 (1985); Phys. Fluids 29, 341 (1986). 

[8] C. T. Hsu, R. D. Hazeltine and P. J. Morrison, Phys. Fluids 29, 1480 (1986). 

[9] A. J. Brizard, Phys. Fluids B 4, 1213 (1992). 
[10] A. B. Mikhailovskii and V. S. Tsypin, Plasma Phys. 13, 785 (1971). 
[11] F. L. Hinton and C. W. Horton, Phys. Fluids 14, 116 (1971). 
[12] B. Scott, Phys. Plasmas 10, 963 (2003). 

[13] R. D. Hazeltine, C. T. Hsu and P. J. Morrison, Phys. Fluids 30, 3204 (1987). 

[14] T. Eickermann and K. H. Spatschek, Phys. Plasmas 3, 2869 (1996). 

[15] O. E. Garcia, J. Plasma Phys. 65, 81 (2001). 

[16] D. Dagnelund and V. P. Pavlenko, Phys. Scr. 71, 293 (2005). 

[17] P. J. Morrison, I. L. Caldas and H. Tasso, Z. Naturfors. Sect. A-J. Phys. Sci. 39, 1023 (1984). 
[18] S. I. Braginskii, Rev. Plasma Phys. 1, 205 (1965). 



